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I-    INTEODOCTION 

There  are  many  solutions  to  the  problem  of  predicting 
the    location     of   an    underwater      energy   source.  One   common 

solution  is  the  use  cf  a  passive  hydrophone  which  detects 
the    pressure  waves      radiating    from  the  source.  The   hydro- 

phone sensor  is  assumed  to  be  omnidirectional  and  therefore 
incapable  cf  estimating  direction.  To  provide  direction- 
ality a  series  cf  sensors  are  placed  in  a  row  to  form  a 
passive    linear    array. 

A  familiar  method  of  determining  dir actionality  is 
time-domain  beam  forming.  In  this  principle,  it  is  assumed 
that  the  source  is  far  enough  away  so  that  the  pressure  wave 
appears  to  be  a  plane  wave  when  viewid  at  the  site  of  the 
receiving   array. 

Thus  a  set  of  time  delays  are  calculated  for  any  direc- 
tion cf  signal  arrival,  which,  when  applied  to  the  receiver 
outputs  causes  them  to  be  in  phase  and  to  reinforce  when 
summed.  -The  resultant  angular  response  to  signals  arriving 
from  other  than  the  nominated  direction  is  then  a  function 
cf  the  array  geometry,  relative  to  the  signal  wavelength, 
and  any  weighting  factors  which  have  been  applied  to  the 
receiver      outputs.  The      effect     is      to      generate      a      main 

receiving  beam  in  the  desired  diraction,  with  a  series  of 
undesired  subsidiary  sidelobes  whose  magnitude  can  be 
controlled  to  some  extent  by  the  choice  of  a  suitable  array 
geometry  and  the  use  cf  amplitude  weightings  on  the  receiver 
outputs. 

In  order  to  determine  location,  three  or  more  such 
arrays  separated   by    a  known    amount   may   be   used. 

This  study  is  concerned  with  a  single  linear  vertical 
passive    array  and  the   determination   of   the    depth   of   a   single 


underwater   source.         The   analysis   is      based  on  the    following 
assumptions: 

•  The  underwater  source  is  emitting  continuously  and  at  a 
monochromatic  frequency. 

•  Both  the  source  and  receiving  hydrophones  are  stationary 
in   space  causing  a  constant   rangs. 

•  The  range  is  sufficiently  long  so  that  the  channel  is 
filled   with    R-R     (ref  rac  ted-refrac  ted)    rays. 

•  There  is  no  distortion  introduced  in  the  propagating 
medium  so  that  the  signals  received  at  each  sensor  are 
identical   except    for  constant  dalays. 

•  The  source  signal  and  noise  are  independent  and 
stationary    gaussian  random    processes. 

•  The  speed  of  sour.d  profile  is  triangular  and  symmetric 
with  the  deep  sound  channel  axis  at  1000  meters.  The 
velocity  gradient  is  -0.017  maters/meter /sec  above  1000 
meters  and  +0.017  meters/meter/sec  below  1000  meters. 
This  profile  gives  a  speed  of  sound  at  ^he  surface  of 
1500    meters/sec. 

•  The  speed  of  sound  profile  is  constant  in  the  horizontal 
plane. 

•  Only  R-R  rays  are  considered.  All  other  rays  have 
sufficient    less    that   their  effect   is   negligible. 

As  opposed  to  conventional  time-domain  beamf orming,  this 
study  makes  no  assumption  of  planar  wave  fronts  at  the 
receiver  site.  Therefore  the  time  delays  applied  to  each 
receiver  will  not,  in  general,  be  a  linear  function  of 
depth. 

Since  it  is  desired  to  determine  whether  or  not  there  is 
a  source      present  at    a  specific      depth   the    result    will      be   a 
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binary      decision.  A     "1"      will     indicate     signal      source 

present.;  a  "0"  will  indicate  signal  source  not  present.  For 
the  constant  range  there  will  be  "N"  test  depths  investi- 
gated for  the  signal  source.  The  number  of  hydrophones  in 
the    vertical  array    will   be    "L". 

For  a  single  source  at  a  given  depth,  the  travel  time  is 
calculated  from  the  depth  to  each  hydrophone.  This  -ravel 
time  is  converted  into  a  phase  delay  for  each  hydrophone  so 
that  after  summation  from  all  hydrophones  a  maximum  output 
is  achieved.  This  output  is  then  passed  through  a  squaring 
device,  an  integrator,  and  a  threshold  and  flip  flop  device 
to  give  a  "1"  binary  output.  If  the  signal  source  is  at  a 
different  depth  and  the  same  previous  phase  delays  are  used 
for  each  hydrophone  then  the  output  will  be  somewhat  less 
than  the  previous  maximum.  The  difference  in  depth  required 
to  achieve  a  "0"  binary  output  is  the  depth  resolution  of 
the    system. 

The  travel  times  for  each  hydrophone  are  calculated  for 
each  cf  the  "N"  source  depths  to  be  considered.  "N"  will 
ordinarily  be  much  greater  than  "L"  so  that  the  system  will 
be  overdetermined.  An  over  det  ermined  system  is  one  in  which 
there  are  more  equations  than  unknowns.  The  objective  then 
is  to  calculate  the  phase  angle  and  amplitude  weight  for 
each  hydrophone  so  that  a  determination  can  be  made  indi- 
cating the  presence  or  absence  of  a  signal  source  at  a  given 
depth. 

The  method  used  to  calculate  the  phase  and  amplitude 
weights  is  the  linear  minimum  variance  estimation  technique, 
linear  minimum  variance  estimators  are  optimum  when  compared 
with  all  ether  estimators  for  gaussian  problems.  The  method 
is   directly   applicable  to   overdetermined   systems. 

The  output  of  the  summer  is  calculated  using  the  linear 
minimum  variance  aiplitude  and  phase  angles  assuming  a 
source   at  one  of  the    "N"   source      depths    and  no  source  at   the 


11 


ethers.  The  calculation  is  repeated  for  each  of  the  depths. 
The  result,  when  plotted  against  source  depth,  will  be 
referred  to  as  the  "beam  pattern"  of  the  array  in  this 
report.  (Although    similar,    it    should   not    be   interpreted   as 

the  angular  response  of  an  array  as  in  the  conventional 
definition  cf  a  beam  pattern.  The  conventional  definition 
loses  much  of  its  utility  when  the  wavef rents  are  not 
planar.)  Ideally  the  beam  pattern  will  be  maximum  at  the 
desired  depth  and  very  small  at  ail  other  depths  so  tha4:  the 
binary  "1"  decision  will  be  made  for  a  source  at  the  desired 
depth,  and  a  "0"  decision  for  sources  at  all  others.  This 
team  pattern  is  compared  with  the  depth  beam  pattern  of  the 
conventional  time-domain  beamformer  mentioned  above.  The 
purpose  cf  this  thesis  is  tc  determine,  as  an  initial  inves- 
tigation, whether  the  linear  minimum  variar.es  estimation 
technique,  when  applied  to  a  linear  vertical  array,  is 
useful  in  depth  discrimination  at  long  ranges  in  a  'SOFAR' 
type    sound   channel. 
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II.    GEHERAL    THEORY 

A.       EAT    ACOUSTICS 

The  propagation  cf  sound  in  an  elastic  medium  can  be 
described  mathematically  by  solutions  of  the  wave  equation 
using  the  appropriate  boundary  and  medium  conditions  for  a 
particular  problem.  The  wave  equation  relating  the  acoustic 
pressure  'p'  to  the  coordinates  •x,,'y,,,z,r  and  the  time 
't1  ,    may   te    written    as 

a2  2     A2  A2  A2 

dp         c      dp    +    ap    +    dp 

7~2  T"2"         7~2  T~2  (2.1) 

dt  dx  dy  dz 

where  'c*  is  a  quantity  that  has  the  general  significance  of 
sound   velocity   and    may  vary    with   the    coordinates. 

One  may  approximate  the  solution  of  the  wave  equation 
using  .ray  theory:  its  body  of  results  and  conclusions  is 
called   ray   acoustics. 

Officer  [Hef.  1]  describes  the  ray  solution  as  a 
complete  solution  tc  any  particular  propagation  problem 
within  the  validity  cf  the  approximation  of  the  Eikonal 
equation  to  the  wave  equation.  For  these  approximations  to 
be  valid  neither  the  amplitude  of  the  wave  nor  the  speed  of 
sound  can  change  appreciably  in  distances  comparable  to  a 
wavelength. 

Thus  the  path  of  a  ray  through  a  medium  in  which  the 
speed  of  sound  varies  with  depth  can  be  calculated  by  the 
application   of    Snell's  law 


cos6/c   =    1/c      =  a    constant   for   any    one   ray  (2.2) 
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where  'e1  is  the  angle  of  depression  made  with  the 
horizontal  at  a  depth  where  the  speed  of  sound  is  '  c*»  and 
•cq '  is  the  speed  at  a  depth  (real  or  extrapolated)  where 
the   ray    wculd  become    horizontal. 

In  a  medium  in  which  the  velocity  of  sound  changes 
linearly  with  depth  the  sound  rays  can  be  shewn  to  be  arcs 
of  circles,  that  is,  to  have  a  constant  radius  of  curvature. 
Kinsler  et  al.  [Ref.  2]  give  a  simple  and  heuristic  demons- 
tation  of  the  circularity  of  rays  in  a  medium  with  a  linear 
sound  speed  gradient  'g'.  The  center  of  the  circle  which 
creates  the  arc  lies  at  a  depth  where  the  sound  speed  extra- 
polates to  zero.  To  understand  this,  consider  a  portion  of 
a  ray  path  with  a  local  radius  of  curvature  'R',  as  illus- 
trated in  Figure  2.1.  Since  the  gradient  'g1  for  this  case 
is 


g=Ac/Az   =  ^2-0^  /  (d2~&1)  =  (€2-0^  /R(cosQ1-cosQ2)  (2.3) 

where  *  a  c*  is  the  change  in  sound  speed  and  'a  z'  is  the 
change  in  depth.  It  can  be  seen  that  the  radius  of  curva- 
ture   is    given  by 


R    =   -c   /g   =  -c/(g    ccse)  (2.4) 

The  ray  path  is  therefore  a  circle  when  'g*  is  constant 
because  *R'  is  then  constant.  The  center  of  curvature  of  a 
circle  lies  at  the  depth  where  6  is  90  degrees,  which  corre- 
sponds to  c=0.  For  the  situation  in  Figure  2.1  the  speed 
gradient  is  negative  so  that  '  R'  is  positive.  If  the  speed 
gradient  were  positive  • R*  would  be  negative,  and  the  path 
wculd   curve    upward. 

Once  the  radius  cf  curvature  of  each  segment  of  a  path 
is      known      the    actual     path      can      be    traced      graphically     or 
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Velocity 

profile 


Figure  2.1   Circular  Hay  Path 
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computed.  If  the  initial  angle  of  depression  of  a  ray  is 
G1 ,  then  by  referring  to  the  geometry  of  Figure  2.1  the 
changes   in    both    range  and  depth   are 

Ar=   c-l  (sin  (9  1)  -sin  (92)  )/(g  cos(91))  (2.5) 

A  z-   c±  (cos  (92)  -cos  (91)  )/(g   cos  (91)  )  (2  .6) 

The  sign  convention  for  these  equations  is:  downward,  to  -he 
right,  and  depression  angles  below  the  horizontal  axis 
positive. 

The  symmetric  triangular  sound  speed  profile  assumed  in 
the  introduction  is  similar  to  speed  profiles  encountered  in 
the  deep  sound  channel,  sometimes  called  the  SOFAR  channel. 
The  velocity  minimum  which  occurs  at  the  axis  of  the  sound 
channel  causes  the  sea  to  act  like  a  kind  of  lens;  above  and 
below^the  minimum,  the  velocity  gradient  continually  bends 
the  scund  rays  toward  the  depth  of  minimum  velocity.  A 
portion  of  the  power  radiated  by  a  source  in  the  deep  sound 
channel  accordingly  remains  within  the  channel  and  encoun- 
ters nc  accustic  losses  by  reflection  from  the  surface  and 
bottcm.  These  rays  are  called  R-R  (r ef racted-ref racted)  and 
since  they  have  very  low  transmission  loss,  very  long  ranges 
can  be  obtained  from  a  source  of  moderate  acoustic  power 
output.  Thus  an  energy  source  at  a  specific  depth  will 
propagate  energy  in  all  directions  but  only  the  direction 
which  is  toward  the  receiving  array  and  for  which  the  rays 
are  R-R  is  of  interest.  To  determine  the  range  of  depres- 
sion angles  which  will  yield  R-R  rays  Snell's  law  is  used  to 
give 
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9raax   =   arc  cos  (c1/  (d-  (g    d) )  )  (2.7) 

where  • emax'  is  in  radians  and  •c1l  is  the  speed  of  sound  at 
the   scurci    depth    ■ d'- 

An  example  is  given  in  Figure  2.2  of  a  single  ray  trace 
propagating  in  the  SOFAR  channel  lo  show  how  equations  2.5 
and  2.6  are  used  in  determining  range  and  depth.  The  ray 
path  is  broken  up  into  arcs  of  circles  as  shown  in  -che 
figure,  and  then  by  paying  close  atten-ion  to  the  previously 
defined  sign  conventions,  the  change  in  range  and  depth  is 
found.  Being  more  specific,  for  arc  1,  91  and  62  are  both 
positive;  for  arc  2,  61  is  positive  and  62  is  zero;  for  arc 
3,  ei  is  zero  and  62  is  negative;  and  for  arc  4,  61  is 
negative  and  82  is  zero.  3y  keeping  a  running  total  of  all 
the  depth  and  range  changes  it  is  possible  to  determine  the 
total  horizontal  distance  travelled  and  the  depth  at  that 
distance. 

For  the  spaed  of  sound  profile  assumed  in  the  introduc- 
tion a  computer  generated  ray  plot  is  shown  in  Figure  2.3 
for  a  source  depth  of  300  meters.  As  each  ray  propagates 
out  from  the  source  the  triangular  channel  becomes  filled 
with  sound.  -If  a  receiving  hydrophone  is  placed  a  great 
distance  away,  a  number  of  refracted  propagation  paths  will 
exist,  each  having  a  different  travel  time  and  crossing  the 
channel  axis  at  different  intervals.  The  path  with  the 
greatest  excursion  from  the  axis  will  have  the  shortest 
travel   time. 

Officer  [Ref.  1]  shows  that  the  travel  time  *  t'  of  a 
ray,    which   is  an  arc   of   a  circle,    is    given   by 


i*1 

g   0     cos9  (2-8) 


t  =  i  c2  de 


i 


and   the   travel   time    for   each   arc   is 


17 


t 


1    log    ftan(  (tt/4)  +  (92/2)) 
=    "    g  e[tan((Tr/4)  +  (61/2))> 


(2.9) 


Equation  2.9,  appliad  using  the  same  convention  as  equations 
2.5  and  2.6,  determines  the  total  travel  time  of  a  ray  in 
the   deep   sound    channel. 

B.       SBHAY    HCDEL 

The  receiving  linear  vertical  array  is  presumed  to 
consist  of  'L'  hydrophones  as  shown  in  Figure  2.4.  It  is 
assumed  that  the  source  is  emitting  energy  at  a  constant 
frequency,  ' f * ,  and  amplitude,  'A',  regardless  of  the  depzh. 
The      source     signal    at     the      source      is    Aexp(j21T    ft).  The 

inherent   received  signal  at    the   firs-    hydrophone    is 


x1(t)    =   A'exp  (  j2^f  (t-t1)) 


(2.  10) 


where  't1'  is  the  travel  time  from  the  energy  source  to  the 
first  hydrophone  and  fA*  ■  is  the  amplitude  of  the  signal  at 
the  range  cf  the  array.  After  passage  through  the  amplitude 
weight  'a1*  and  a  phase  delay  of  '  ti«,  the  signal  on  the 
first   hydrophone   at    the   input   to   the    summer   is 


yl(t)=A'alxl(t-Ti  )=A'alexp(  j2-rrf  (t-tl-n  )) 


(2.11) 


A  ^ime   delay,    for  monochromatic   signals   corresponds  to  a 
phase    shift 


9    =    2^  ft 


(2.  12) 


where    •9'    is  the   phase  shift   in      radians   and    't»    is    the   -ime 
delay   in   seconds.      Thus    equation    2.11    can    be    written   as 
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f  1  (t)  =A»a1sxp(  j(2    ft-<j>1-91))  (2.13) 

where  '  §  1=2irft  1  *  is  the  phase  delay  due  to  the  travel  time 
from  the  source  re  the  first  hydrophone  and  '91'  is  the 
phase   dalay   in    the    receiver    on    the   first   hydrophone. 

Combining  all  the  hydrophones  in  the  array  in  a  sumier 
gives   as   an   expression  for    the   array    output 

L 
Y(t)    =     I    &«akaxp(j<2irft-fk-9k))  (2.14) 

where  '  <j>k  •  represents  the  phase  delay  due  to  the  -ravel  time 
from  tha  source  to  the  "k  th"  hydrophone  and  *  9,  •  is  the 
phase  delay  in  the  receiver  on  the  "k  th"  hydrophone.  If 
the  amplitude  of  the  energy  sourca  is  normalized  by  setting 
A'  =  1  ,  and  equation  2.14  is  written  in  terms  of  real  and 
imaginary   components,    we   have 

L 
Y(t)  =      5>k   cos(-(J>k-8k)+jsin(-4>k-ek)    exp(j2Trft) 

When  an  energy  source  is  at  the  "q  th"  depth,  we  wish  to 
have  each  receiving  hydrophone's  phase  delay  cancel  out  the 
effect  of  the  travel  time  from  tha  source  to  it,  such  that 
'"fk'^k1  ^s  squal  to  a  multipla  of  ,2^' .  This  will  put  all 
signals  into  the  suirmer  in  phase  and  thus  maximiza  the 
signal  gain  for  a  source  at  the  "q  th "  depth.  Frcm  equation 
2.15, 

L 

I    (a,   cos("(^)]ca~eka')  ^  =  L   ^source  present  at  q)   (2.16) 


I    (a.   sin(-<£   -8   ))  =  0   (source  present  at  q) 
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Nots  that  the  firs-  subscript  on  -he  phase  angle  indicates 
the  receiving  hydrophone  and  the  second  subscript  indicates 
the  depth  of  the  source.  Thus  ' 6^  '  would  indicate  the 
phase  shift  relating  to  the  travel  time  from  the  "q  th"  test 
depth   to   the   "k    th"    hydrophone   in   the   receiving   array. 

It  is  desirable  for  'Y(t)  •  to  be  a  minimum  value  for 
sources  at  other  than  the  depth  being  investigated.  Thus 
for  each  of  the  other  *N-1'  depths  'Y  (t)  *  is  set  to  zero. 
This  gives  'N-1*  equations  for  the  real  terms  of  •Y(t)  '  set 
to    ze  ro 

y    a^qcos  (-^km'^a)  =0    (source    absent",    m/q)  (2.18) 

k=l 

and  'N-1'  equations  for  the  imaginary  terms  of  'Y(t)  '  set  to 
zero 


L 

k 


I    ak   sin  (-f^-^g)  =°      (source   absent;    m/q)  (2.19) 


By  using  elementary  trigoncmetric  identities,  equation 
2.16  (real  terms  with  source  present  at  "q  th"  depth) 
tecomes 


L 
k 


I       akq  fcos  (fkq)cos  (Qkq)-sin(^q)5in  (Okq)j  =L  (2.20) 


Equation    2.18    (real    terms   with  source   absent    for    each   of   the 
ether    'N-11    depths)     becomes 

I   akq[cos(*km)cos(9kq)-sin(*km)sin(9kq)]=0   m=l,2,...N;    m/q  (2^21) 

Equation    2.17    (imaginary   terms   with   source   present   at   the   "q 
th"    depth)    becomes 
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y   -a,  (sin  (6      )cos(6     )+ccs(4      )  sin  (9     )  )  =0 


k=l 


kq 


kq  kq' 


kq  kq 


and   equation      2.19    (imaginary   terms      with   source 
each    of   the   other    •N-11    depths)     beccmas 


(2.22) 

absent    fo: 


L  - 

J1-akq[sin(*km)cos(9kq>+c°sC*kni)sin(ekq)]  =0 


i=l,2,  ...N;    m^< 


(2.23) 


Thus,        there      are      a     total   of      '2N'       equations      with      '  2L' 
unknowns. 

In  order  to  simplify,  we  put  these  '2N'  equations  into 
matrix  fcrm.  Arbitrarily  the  real  terms  are  made  the  first 
•N»  equations  and  the  imaginary  tsrms  the  sscond  '  N'  equa- 
tions. The  first  real  and  first  imaginary  equation  is  at 
the  lowest  (shallowest)  source  depth  and  equations  increase 
in  order  after  that  until  the  last  real  and  last  imaginary 
equation  correspond  tc  a  source  at  the  deepest  depth.  The 
resultant   matrix  equation   becomes: 


cos4>11 

cos4>21 

•    cos<f,Ll 

-sin*,. 

-sin<t>21    .  . 

•    -sin*Ll 

cos<(>12 

cos<*22 

•    COS(fL2 

-sin4>12 

-sin4>22    •  • 

•    -sin*L2 

cos*lq 

cos02q 

•    cos*Lq 

"sin*lq 

-sin<t»2q    .. 

•    -sin0Lq 

eos$1N 

GOS(j)2N 

•    GOS,*LN 

-sind.1N 

~sin*2N    *" 

•     -sin*LN 

-sin<fr11 

-sin* 
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.  -sin*.  , 

-cos<()11 

-cos4>21    . 

.    -cos*. . 

-sin<J>12 

-sin* 
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.-sin*L2 

-cos*12 

-cos<i)22     • 

•    -cos*L2 

-sin*1N 

-sin* 

2N 

.-sin*LN 

-cos<(.1N 

-cos*2N    . 

'    -cos*LN 

alqCOs9lq 
a2qcos92q 


a.  cosQ. 
Lq     Lq 


alqsin9lq 
a2qsin92q 


a.  sin 


[_Lq   LqJ  (2.  24) 


Simplifying    further,    equation    2.24   becomes 
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L 


0 


[       " 

a 

-8 

-3 

-a 

_ 

~4 

(2.25) 


where    •    a.1,      •_§_'  ,      'cj ,        and    'sj    represent      the    appropriate 
submatrices. 

Then,  ty  noting  than  the  multipl icaticn  of  each  element 
of  a  column  by  the  same  nonzero  constant  doesn't  affect  the 
solution,   equation    2.25   becomes: 


0 


a 

e 

-i 

a 

0 


Finally   letting    Z,    A,     and    0   represent 
tion    2.26,    we  obtain 


(2. 26) 


the    matrices    in    equa- 


Z    =    A0 


(2.  27) 


where      a      matrix     is    denoted      ty      a      capitalized      underlined 
letter. 

In  summary,  Z  is  the  desired  response  of  the  vertical 
array  to  the  'N'  source  test  depths.  A  represents  known 
travel  times  from  each  of  the  'N'  source  depths  to  each  of 
the    '  L'    receiving      hydrophones.       J2.  is   unknown.  It   is   the 

phase   and  amplitude      weighting   which    must    be     applied   tc   the 
vertical     array    in      crder  to     realize   T±.  8   contains      '  2L' 

unknowns. 
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Equation  2.27  represents  '2N*  equations.  Since  this 
system  of  equations  is  over  determined  (N>L)  an  exact  solu- 
tion dees  net  exist.  In  order  to  make  the  best  estimate  of 
6  for  the  .desired  response,  the  linear  minimum  variance 
estimation    technique    is   used. 

C.       LINEAB    HINIMDM    71BIANCE    HETHOD 

Equation  2.2  7  represents  a  noise  free  environment.  If 
noise    ware    present    it   would    become 

_Z    =    A9+n_  (2.  28) 

with  ,n>1  a  "2N"  element  column  vector.  This  represents  the 
noise  at  each  source  depth.  '_Z_'  is  a  linear  function  of  '6' 
and   is   called  the  oberserva ticn.  *kj    is   a    "2Nx2L"   modula- 

tion or  observation  matrix  which  is  known,  and  '.e*  is  the 
"2L"  element  random  parameter  vector  which  is  to  be  esti- 
mated. Assume  the  first  and  second  moments  of  'JS*  and  'n/ 
are    qiven   by 

E(iL)    =  iie  Var(9)    ■    VQ  (2.29) 

and 

E(n)    =0  Var(n)     =   Vn  (2.30) 

where  'E*  represents  expectation  or  first  moment  and  'Var' 
represents  the  covariance  matrix.  It  is  assumed  that  the 
parameter    'G1   and  the  noise    'n.'    are   uncorrelated. 

A  restriction  imposed  is  that  the  estimate  must  be  a 
weiqhted   linear    combination    of  the  observations: 
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§T     ■    b+£Z.  (2.31) 

where  ' * ■  indicates  estimate.  The  objective  is  re  select 
'bf  and  »JB/  in  order  tc  minimize  the  error  variance.  Such 
an  estimator  is  the  linear  minimum  variance  estimator;  it  is 
the  best,  in  the  sense  cf  minimum-error-variance  linear 
estimators. 

Another  restriction  is  that  the  estimator  be  unbiased; 
in  otherwords  it  is  required  that  the  expected  value  of  the 
estimator    »e  •    is  equal   to      the   expected    value  of    the   param- 

Li 

eter    '0'.      Thus, 


E(eL)=b  +  BE(z)=E(Q)  =    yQ  (2.32) 

yielding 

b=yQ-EAUQ  (2.33) 

Substituting      this    result      in    equation      2.31      gives    for      the 
unbiased   linear    estinator 


©L=yQ+E(Z-AyQ)  (2.  34) 

Note  that  since  the  estimator  is  unbiased,  the  estimation 
error  '9  =6-6  ■  is  zero  mean.  The  next  step  is  -c  select 
'_3'  in  order  to  minimize  the  error  variance.  However,  this 
optimization  problem  is  ill-defined  because  the  error  vari- 
ance is  a  matrix.  Therefore  in  order  to  introduce  a  scalar 
goodness  measure  the  sum  of  the  variances  of  each  component 
of  'JM  is  minimized.  This  is  the  sum  of  the  main  diagonal 
terms  of  the  co variance  matrix  and  is  defined  as  the  trace 
cf    the   matrix. 
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2N 

tr    (Var(9r))    =         Y      Var((9_.)    )  (2.35) 

— E  u  n  — E    n 

n=l 

where      *tr'      indicates     trace.         'B*      is      than      selected     to 
minimize   the  trace    of  the  errcr  variance,    cr 


min    tr    (Var(0    ))     =    min   tr     (E  (6   eS  )  (2.36) 

B  ~L  B  L-  L 

T 

where  ■ x  '  indicates  the  transpose.  The  following  problem  is 
then  obtained  by  substituting  equation  2.34  into  equation 
2.  36  : 


min   tr    ( Var(9„) )  =min   tr(E(  (9-u_ -3(Z-AjO  (9-ua3(Z-Au~  )T) )        (2.37) 

_  —  c  _  —  —o    —  —  — — o       —  — o  —  —  d  *  ' 

It    is      well    known  [  Hef •    3]      that  equation    2.37      is  minimized 
when 


Cov  (9,Z)-BVar  (Z)  =0  (2.38) 

where    '  Ccv(e,Z)  '      is    the      ccvariance  matrix      of   the      unknown 

parameters      and      the   observations.  Denoting     the      optimum 
filter    by    f.B:':?  \    then   if 

_B*=Cov  (6,  Z)  (Var(ZJ  )  _1  (2.39) 

a  minimum      is  achieved      for   the      sum  of      the   squares      of  the 
errors. 

Using   equation    2.28   for    %Z} ,        the   co variance    of    '.9J    and 
%Z}    becomes 

Cov  (0,Z)=Cov  (6rAG  +  n)=V6AT  (2.U0) 

since    •e.'    and    'jo,1    are  un correlated.  The   variance    of    '_ZJ    is 
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Var  (ZJ  =Var  (A9+nJ  =£V   AT+V  (2.  41) 

Substituting  these    into   equation  2.39    gives 

5*    =    Z8^T(^-9^T+-n)"1  (2.42) 

and   the    linear    minimum  variance   estimator   is 


iLMV    '   tVIe^AVeA^r^Z-A^)  (2.U3) 

By  utilizing  a  matrix   inversion  lamma  [Ref.  3]  equation 
2.43  b€ccmes 


iLMV  "  (ATi;1A+^1)_1(ATv;1Z+V^1!i8)  (2.  M) 

The  advantaqe  of  this  equation  over  equation  2.43  is  the 
size  of  the  matrix  to  be  inverted.  In  equation  2.4  3  the 
matrix  has  dimensionality  • 2N *  whila  in  equation  2.44  its 
dimenionality  is      only   • 2L1  .  Thus    the      advantages   of      the 

linear  variance  estimator  are  the  ease  with  which  they  are 
derived,  the  mathematical  tractability  of  the  linear  form, 
and  the  minimum  amount  of  stochastic  information  required 
for  development.  An  interesting  characteristic  is  that  the 
linear  minimum  variance  estimate  is  the  orthogonal  projec- 
tion of  '9*  onto  the  space  spanned  by  the  observation  ' Zf . 
Eecause  of  these  factors  this  estimator  is  a  popular  form 
for    estimating   unknowns   in    cverdetermined   equations. 

For   this  thesis    it  is   assumed   that   the   noise    samples   are 
uncorrelated  and  identically   distributed  so   that: 


V    =  gl  (2.45) 

~~n 
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No    previous   knowledge   is   assumed   about    '9'.      This    implies  an 
infinite    variance   matrix   which   is   represented   as: 


7n    =0 


and 


U0  =0 


(2.46) 


The    linear    minimum    variance    estimate      given   by  equation    2.44 
is   then 


£utf<iT«"Vi 


(2.47) 


By   determining    •  G      ■    the  phase   and   amplitude    weights   are 
found   for   a   signal    source   on  the    'q   th'    depth.      Recall   that 


-LMV 


'L+l 


alqcos9lq 

a2qCOs92q 


aLQCOS9Lq 


alqSlnaiq 


a.  sine. 

-q   ^q 


(2.48) 


2L.    L 

Dpon    solving,   this    equation    gives   for   the   phase   delay 


9      =arctan(9T    /9   )    where    m=1   to   L 

ma  L+m  m 


(2.49) 


and    amplitude   weight 


a      -9   /cos  (9      )  where    m=1    to   L 

mq      m  mq 


(2.50) 


When  these  amplitude  weights  and  phasa  delays  are  applied  to 
the  vertical  linear  array  a  resulting  beam  pattern  is  formed 
which   in   the  absence    of  noise    is: 
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=  ae 


(2. 51) 


The  resulting  beam  pattern  *Z_  '  ca.n  then  be  compared  with 
the  desired  beam  pattern  ,Z_t  ,  as  wall  as  with  a  conventional 
beam  pattern  'ZJ.1  obtained  using  linaar  phase  shifts  across 
the  array  apertura  selected  to  "steer"  the  array  to  the 
dominant   arrival   angle  for    the   selected   source  depth. 
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III.    EXPERIMENTAL    PROCEDURE 

A.  EASIC    ASSUHPTIONS 

The  speed  of  sound  profile  given  in  the  Introduction  was 
used  to  test  the  techinique ,  with  'L',  the  number  of  source 
depths  chosen  as  20.  The  source  depths  are  at  220  i^ers 
and  every  20  meters  thereafter  to  6  00  meters  inclusive. 
•N1,  the  number  of  hydrophones  used  in  ths  vertical  array 
was  chosen  as  5.  The  hydrophones  are  placed  at  depths  of 
100,  200,  300,  400,  and  500  meters.  Thus,  there  are  40 
equations,  two  for  each  source  depth,  and  10  unknowns,  two 
for  each  hydrophone.  A  range  of  200  kilometers  (km)  assures 
that  the  deep  scund  channel  is  filled  with  sound  over  the 
aperture  of  the  vertical  linear  array.  A  frequency  of  100 
Hz  provides  good  resolution  in  the  beam  pattern  without 
introducing  alias  mainlobes  at  the  selected  range  across  the 
depths   of   investigation. 

B.  «A/     HATRH    CALCULATION 

Since  the  gradient  of  the  sound  velocity,  ,g«  is 
constant  in  the  area  of  the  source  depths,  and  the  speed  of 
sound  at  the  surface  is  known,  equation  2.3  is  used  to  solve 
for  'C21,  the  speed  of  sound  at  the  source  depth.  Equation 
2.7  is  used  for  each  source  depth  to  calculate  the  maximum 
initial  depression  angle  which  yields  R-R  rays.  In  deter- 
mining travel  time,  only  depression  angles  from  each  source 
which   are   positive    (downward)    and   yield    R-R   rays    are  used. 

Beginning  with  the  first  source  depth  of  220  meters,  an 
initial  depression  angle  of  0.0  degrees  is  selected.  The 
ray  path  is  then  calculated  using  aquations  2.5  and  2.6  and 
broken   into   a  series   of   arcs   as   in   figure   2.2.      Equation   2.9 
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is  used  to  determine  the  travel  time  for  each  arc  ir.  the 
same  manner.  Each  arc's  horizontal  range,  travel  time,  and 
depth  are  summed.  When  the  summed  horizontal  range  reaches 
200    km      the   summation     process   ends.  The   travel      time   and 

depth  of  the  ray  at  this  horizontal  range  is  then  known. 
The  same  procedure  is  repeated  for  an  initial  depression 
angle  of  0.1  degrees  and  every  increment  of  0.1  degrees 
thereafter  until  the  naximu m  depression  from  equation  2.7  is 
reached.  The  final  ray  path  is  at  this  maximum  depression 
angle. 

The  same  procedure  is  repeated  for  each  of  the  other  19 
source   depths. 

Thus,  for  each  initial  angle  from  each  source  d^pth 
there  is  a  ray  which  has  a  travel  time  and  a  depth  when  it 
reaches  the  horizontal  range  of  200  km.  Since  it  is  the 
profile  of  the  sound  pressure  wave  which  impinges  on  the 
vertical  array  which  is  of  importance,  a  constant  can  be 
subtracted      from      these     calculated        travel      times.  This 

constant  is  selected  to  be  the  travel  time  for  the  source 
depth  of  220  meters  which  has  an  initial  depression  angle  of 
0  degrees.  It   is    subtracted   from   each  of  the  travel   times 

making  the  resultant  travel  times  relative  with  respect  to 
the  ray  which  has  a  0  degree  depression  angle  from  the  220 
meter  depth.  The  program  and  its  listing  which  calculates 
the  relative  travel  times  and  the  depths  of  these  rays  at 
the    horizontal    range   of   200    km    is   given   in   Appendix   A. 

A  plot  of  the  relative  travel  times  versus  depth  for  the 
220  meter  source  depth  is  shown  in  Figure  3.1.  Figure  3.2 
displays  the  plot  for  the  380  meter  source  depth.  Negative 
relative  travel  times  in  the  plots  indicate  that  the  overall 
travel  time  is  less  than  the  reference.  These  rays  arrive 
at   the    200    km  horizontal  distance   before   the   reference   ray. 

Since  the  receiving  hydrophones  are  at  set  vertical 
positions       (100,      200,        300,      400,         and      500  meters),        an 
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interpolation  is  dons  to  determine  relative  travel  times  to 
them  from  each  source  depth  .  The  interpolation  program  and 
its   listing   is    in   Appendix    B.  Sometimes    mors   than   one   R-R 

ray  travels  from  the  source  depth  to  a  hydrophone.  When 
this  cccurs,  the  ray  which  arrives  first  is  used  in  the 
calculation   of    relative  travel  time   to   that   hydrophone. 

Eguation  2.12  determines  the  phase  shift  relating  to  the 
relative  travel  times.  The  *>AJ  matrix  is  formed  by  taking 
the  appropriate  sine  and  cosine  values  as  in  equation  2.24. 
The    ,JP    matrix    is    «40   by    101. 

C.  •£■    HATBIX 

Referring  to  equation  2.26,  the  *ZJ  matrix  is  a  *  40  by 
1'  column  vector.  It  is  the  desired  beam  pattern.  The 
bottom  20  rows  give  the  imaginary  tarms  and  are  set  to  zero. 
The  top  20  rows  represent  the  value  of  the  real  terms  at 
each  source  depth.  Therefore  each  of  the  top  20  rows  is  set 
■to  zero  except  for  the  row  containing  the  source.  It  is  set 
to  1.  Fcr  example,  if  the  source  is  at  220  meters  then  only 
the   tcp   rcw     is    set    tc   1 .  If   ths   source   is      at    380   meters 

then    only  the  ninth   row   is    set   to    1. 

D.  BES0TTIBG    BE1H    PITTEBN 

1  •      S.§ilL2  the  Linear    Minimum   Var i  ance   Met  hod 

1 6TW„   ■     is      calculated      using      equation    2.47.  The 

— LMv 

resulting  beam  pattern  'Z^1  is  calculated  using  equation 
2.51.  The  program  which  calculates  the  'A'  mar.rix,  uses  it 
in  determining  'g,^',  and  then  calculates  'Z_'  is  given  in 
Appendix   C.      The   program   listing   is   also   included. 
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2 •      2§iHa  Linear    Phase    Shifts 

The  conventional  beam  pattern  is  determined  by  using 
equation  2.51  where  '©.'  is  calculated  by  approximating  the 
plot  of  relative  travel  time  vs.  depth  by  a  straight  line  at 
the   receiving     hydrophone  depths.  For   example      Figure   3.3 

represents  this  plot  for  the  380  meter  source  depth.  The 
straight  line  is  determined  by  a  least  squares  linear 
regression  which  minimizes  the  sum  of  the  squares  of  the 
deviations  of  the  actual  data  points  from  the  straight  line 
of   best    fit.  Note   that  only  data    points    which      are   on   the 

dominant  curve  are  used  in  calculating  the  straight  lire. 
From  the  straight  line,  relative  travel  times  to  the 
receiving   hydrophones  are  calculated.  The   relative   travel 

times  for  the  3  80  meter  source  depth  are  given  in  Table  I. 
They  correspond  to  a  plane  wave  arrival  angle  of  3.73 
degrees.  •9J  is  determined  by  converting  these  relative 
travel  times  to  phase  delays  using  equation  2.12  and  then 
taking  the  appropriate  sine  and  cosine  values  of  these  phase 
delays      as   in     equation    2.24.  The      amplitude   weights      are 

initially   assumed  to   be   unity. 

A  second   method   for      obtaining  the   conventional  beam 

pattern      is   calculated     by      the      same    procedure      except      the 

amplitude    weights   which   are    determined   by    equation   2.50,    the 

•9         '      amplitude  amplitude      weights,      are      applied   to      each 
~LMV 
hydrophone. 
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TABLE    I 

Relative   Travel  Times    For   380   Heter   Source    Depth 


Hydrophone    Eepth 

100  meters 

200  meters 

500  meters 

U00  meters 

500  meters 


Relaoivs  Travel  Time 
-0.07565509  sec. 
-0.07131599  sec. 
-0.06697690  sec. 
-0.06263780  sec. 
-0.05829871    sec. 


i 
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IV.     RESOLTS 

A.  EXACT    SOIOTION 

An  sxac-c  solution  is  derived  for  the  beam  pattern  if  the 
number  of  receiving  hydrophones  equal  the  number  of  source 
depths.  For  example,  when  the  5  receivers  are  used  to 
discriminate  between  5  source  depths  (220,  240,  260,  280, 
and  300  meters)  there  are  10  equations  with  10  unknowns. 
Figure  4.1  is  a  plot  of  the  resulting  beam  pattern  with  the 
source  at  the  shallowest  depth.  Note  that  because  of  roun- 
doff errors  in  the  •IMSL'  subroutines  there  is  a  small  value 
for  the  resulting  beam  pattern  at  the  non -energy  source 
depths   under  investigation. 

B.  FOUR    DEPTHS    WITH    TWO    RECEIVERS 

Fcr   source    depths  at    220,  280,    300,      and    320    meters   with 

receiving   hydrophones     at   100  and      200   meters,      there      are   8 

equations   with      4   unknowns.  Figure   4 . 2    is     a    plot      of   the 

resulting   team    pattern   with    a  source    at    the   220   meter   depth. 

Figure      4.3   is      the    plot      for  the      source   at      the    280      meter 
depth. 

C.  TWENTI  DEPTHS  WITH  TWO  RECEIVERS 

For  all  20  source  depths  with  receiving  hydrophones  at 
100  and  200  meters,  there  are  40  equations  with  4  unknowns. 
Figure  4.4  is  a  plot  of  the  resulting  beam  pattern  for  a 
source  at  the  220  meter  depth.  Figure  4.5  is  a  plot  for  the 
source   at   the  380   meter   depth. 

In  order  to  determine  if  the  '91  calculated  in  this  case 
is   the      best  an    alternate      method   is   devised.  Four   source 
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depths    (220,   240,    260   meters,    and   another   source    test   depth) 

with      the      original      source      a-      the      220      meter      depth      and 

receivers   at   100   and    200   meters  are   used.        The   beam    pattern 
v 

is  calculated  each  time  with  a  different  source  test  depth 
substituted     for      the     fourth   source      depth.  The      5      best 

resulting  beam  patterns  are  selected  along  with  the  8  source 
depths  (220,  240,  260  meters,  and  the  5  test  depths  which 
created  the  5  best  team  patterns)  .  Then,  using  these  8 
source  depths,  '9.1  is  determined  for  the  two  receivers. 
This  *QJ  is  applied  to  the  two  receivers  and  the  beam 
pattern   obtained  for   all   20    source   depths. 

The  resulting  beam  pattern  obtained  by  this  alternate 
methcd  isn't  as  good  as  the  beam  pattern  obtained  by  using 
all    20   source  depths   in  the    determination   of    '9'. 

D-       TWEHTY.    DEPTHS    IITB   FITE    RECEI¥ERS 

Fcr  all  20  source  depths  with  all  5  receiving  hydro- 
phones there  are  40  equations  with  10  unknowns.  Figure  4.6 
is  a  plct  cf  the  resulting  beam  pattern  with  the  source  at 
the  220  meter  depth.  Figures  4.7,  4.8,  and  4.9  are  the 
plots  for  the  source  at  the  360,  380,  and  400  meter  depths 
respectively. 

E.       CCHVIBTIOHAL    BEAKFORMER 

Figure  4.10  is  a  plot  of  the  beam  pattern  for  a  conven- 
tional beamformer  using  linear  phase  shifts  across  the  array 
with  the  scurce  at  the  380  meter  depth  and  the  amplitude 
weights  set  to  unity.  All  20  source  depths  and  5  receiving 
hydrophones  are  used.  Note  that  in  figure  4.10  that  there 
is  less  than  1  db  discrimation  between  each  of  the  source 
depths.  Figure  4.11  is  the  plot  obtained  for  the  amplitude 
weights   set   to    values  determined   by   equation    2.50. 
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F.       BAHGE   OF   250    KILCBETERS 

The  calculations  were  repeated  for  a  range  of  2  50  km 
using  the  same  20  source  depths  and  5  receiving  hydrophcnes. 
Figures  4.12  and  4.13  represent  the  plots  cf  relative  travel 
times  versus  depth  fcr  the  220  and  380  meter  source  depths 
respectively.  Figure   4.14      represents      the   straight      line 

approximation  of  the  relative  travel  times  for  the  380  meter 
depth.  Note  that  in  this  figure  the  relative  travel  times 
are  represented  by  two  straight  lines;  the  upper  line  repre- 
sents linear  phase  shifts  of  the  slower  travel  times  for  the 
ccnventicnal  beamforrcer  while  the  low=r  line  represents 
linear  phase  shifts  of  the  faster  travel  times.  Tables  II 
and  III  are  the  straight  line  interpolations  of  these  slower 
and  faster  travel  times  which  correspond  to  arrival  angles 
of  4.04  and  -3.16  degrees  respectively.  Figures  4.15  and 
4.16  are  the  resulting  beam  patterns  for  the  conventional 
beamfcrmer  for  the  slcwer  travel  times  using  unity  amplitude 
weights  and  linear  minimum  variance  amplitude  weights 
respectively.  Figures  4.17  and  4.18  are  the  beam  patterns 
for    the    faster    travel  times. 

Figures  4.19,  4.20,  4.21,  and  4.22  represent  plots  cf 
the  beam  pattern  for  the  source  at  the  220,  340,  360,  and 
380    meter   depths  respectively. 
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TABLE    II 

Slower   Hay    Travel  Ti«es  For   250    Jem   Range  and    380   ■ 

Source 


Hydrophone    Eepth 

100  leters 

200  meters 

300  meters 

400  meters 

500  meters 


Relative  Travel  Time 
-0.09949109  sec. 
-0.09478615  sec. 
-0.09008121  sec. 
-0.08537627  sec. 
-0.08067133    sec. 


■     ■  '   ■-                     ■ ■ 

■* 

TABLE 

III 

Faster  Travel  Tiaes  For    250 

lea 

Range   and  380   m    Source 

Hydrophone    Depth 

Relative   Travel    Time 

100  meters 

-0.  10100745    sec. 

200  meters 

-0.  10468763    sec. 

300  meters 
■ 

-0.  10836781    sec. 

|                    400  meters 

-0.  11204799    sec. 

I                    500  meters 

1 

-0.  11572817    sec. 
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V.    DISCUSS  ION   OF    RESULTS    AND    RECOMMENDATIONS 

For  many  of  the  cases  examined,  the  linear  minimum  vari- 
ance estimation  technique  gives  high  resolution  in  the 
"depth  beam  pattern"  for  sources  at  long  range.  Figures  4. 2 
and  4.3  represent  beam  patterns  for  four  depths  with  two 
receivers  where  the  main  lobe  is  at  the  desired  source  depth 
with      a   beamwidth      of     30   meters.  For      20      depths    with      2 

receivers  the  system  is  highly  overdatermined  and  the  beam 
pattern  portrayed  in  figure  4.4  has  its  main  lobe  80  meters 
from  the  source  depth.  At  the  source  depth  "he  strength  of 
the  beam  pattern  is  3.6  db  down.  For  figure  4.5  the  main 
lobe    is  on   the    source   depth    with   a   beamwidth   of   250    meters. 

For  all  20  source  depths  with  5  hydrophone  receivers  the 
results  range  from  a  beam  pattern  having  its  main  lobe  on 
the  source  depth  (figure  4.7)  with  a  beamwidth  of  20  meters 
and  a  secondary  lobe  12.1  db  down  to  a  beam  pattern  having 
its  main  lobe  40  meters  from  the  source  depth  (figure  4.22). 
For  the  team  pattern  in  figure  4.22  the  width  of  the  main 
lobe  is  80  meters  and  the  strength  of  the  beam  pattern  at 
the   scurce   depth  is    0.8  db    down. 

When  using  the  conventional  beamformer,  the  beam  pattern 
results  with  amplitude  weights  determined  by  the  linear 
minimum  variance  estimation  method  were  superior  to  the  beam 
pattern  determined  by  using  unity  amplitude  weights. 
However  in  all  cases  the  conventional  beamformer  was  signif- 
icantly inferior  to  the  beam  pattern  determined  by  the 
linear  minimum  variance  estimation  technique  in  terms  of 
depth   resolution. 

The  linear  minimum  variance  estimation  technique  could 
not  be  used  for  all  ranges  with  the  chosen  sound  speed 
profile.        For      ranges  of      187    km     and   235      km  the      relative 
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travel   times      produce     an      'A_'      matrix        which     is      so      ill- 
conditioned     that     the      •IMSL'      subroutine      'LINV2F'      cannot 

T 

determine   an  accurate   inverse    of    *k   &*• 

There  is  sufficient  evidence  from  this  initial  investi- 
gation that  the  linear  minimum  variance  estimation  technique 
applied  to  a  long  linear  vertical  array  can  yield  high  reso- 
lution depth  information  about  passive  sources  at  very  long 
ranges.  However,  further  investigation  is  needed  before  any 
field      tests  are      in     order.  Recommendations    for      further 

research   are: 

1.  The  use  of  a  more  realistic  speed  of  scund  profile, 
preferably  one  actually  characteristic  of  the  'SCFAR1 
channel. 

2.  The  use  of  mere  hydrophones  in  the  vertical  array. 
As  more  receivers  are  used  the  beam  pattern  should 
approach  the  desired  beam  pattern  more  closely  (a 
less  over  determined  system  should  yield  smaller  total 
nirimum    mean    square    error). 

3.  Investigation  into  causes  for  the  peak  response 
falling  at  other  than  the  desired  depth  and  techni- 
ques for  correcting  this  problem  when  using  the 
linear  minimum  variance   estimation   techinique. 

4.  Alternate  assumption  for  choosing  the  ray  paths  to 
include  in  the  '_AJ  matrix.  For  example  one  might 
choose  the  rays  with  the  greatest  intensity  instead 
of  the  shortest  travel   time. 

5.  The  possibility  of  estimating  range  as  well  as  -he 
depth  of  a  passive  source  with  linear  minimum  vari- 
ance  estimation  techniques. 

In  conclusion,  the  linear  minimum  variance  estimation 
technique  of  beamforming  was  significantly  superior  to  the 
conventional  beamformer.  High  resolution  depth  information 
about    passive  sources  at   long    ranges   is   provided. 
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APPENDIX    A 
BELATI¥E   TRAVEL    TIHE    CALCULATION 

This  program  calculates  the  relative  travel  times  from 
each  source  depth  tc  the  horizontal  range  and  the  depth  of 
the  ray  at  this  range.  To  save  paper,  only  the  two  shal- 
lowest  source   depths   are   listed. 
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APPENDIX    B 
IHTEBPOLATICN    OF    BELATI7E    TEA7EL    TIME    CALCOLATIONS 

This    program      does  an  interpolation      of   the      output   da*a 
generated   in     Appendix   A.  The   relative      travel    times      are 

interpolated  for  the    receiving   hydrophone    depths. 
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APPENDIX    C 
RESULTING    BEAM    PATTERN    FOR    CALCULATED    WEIGHTS 

This  program  uses  the  output  of  Appendix  B  to  calculate 
the  'A'  matrix.  From  this,  the  amplitude  and  phase  weights 
are  determined  by  the  linear  minumum  variance  estimation 
technique.  These  weights  are  applied  to  the  array  and  the 
beam    pattern  is    obtained. 
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